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A frame- independent, thermodynamically exact London equation is presented, which is especially 
valid for rotating superconductors. A direct result is the unexpectedly high accuracy (~ 10"^") for 
the usual expression of the London moment. 
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The defining property of superconductors to expel any 
magnetic field from its bulk reverses when it is rotated, 
and a spontaneous field (or London moment) 



B 



-2f2/7, 7 = e/mc, 



(1) 



appears, where O is the rotational velocity, while m, e(< 
0) are respectively the bare mass and charge of the elec- 
tron (In MKSA units, 7 = e/m.) This is intriguing 
because some exact cancellations must have taken place 
before microscopic parameters such as e or m will appear 
on a macroscopic level, in a complicated many body sys- 
tem that any superconductor is. So in this sense, the 
London moment may perhaps be compared to the quan- 
tum Hall effect The bare charge e holds its place 
because the superfluid velocity. 



{h/2m)Vip ~ 7A, 



(2) 



must be gauge- invariant. (A is the vector potential, and 
Lp the macroscopic phase variable characterizing the spon- 
taneously broken symmetry.) Although this type of ar- 
guments does not apply to m, a common prefactor in 
Eq(^, many authors have convincingly argued why nev- 
ertheless, and to what extent, m remains unrenormalized 
in the London moment 1^,^-^. 

In this paper, a symmetry argument as simple as the 
gauge invariance is invoked to show why it is the bare 
mass that enters Eqs(|l|,||). The basic point is, Vg is only as 
defined, with the bare mass m, a velocity under Galilean 
transformations. 

There is a preferred coordinate system in superconduc- 
tors, that of the lattice, and most useful properties are 
best understood in this system. Confining one's inter- 
est to these, the transformation property of Vs is irrele- 
vant, and de Gennes' famously pointed statement, that 
the mass m is completely arbitrary, and might well be 
taken as that of the sun [0, valid. 

Sometimes, however, other coordinate systems are 
more advantageous, especially if the lattice itself (when 
rotating or subject to a sound field) is not an inertial sys- 
tem. A theory that consistently relates the results from 
different frames is then desirable. And how Vs transforms 
is an important input for such a frame-independent the- 
ory. As will be shown below, the London moment may be 
understood as a direct consequence of Vs being a velocity. 



The superfluid velocity Vs is a concept perhaps more 
of superfluid helium, and less universally employed in 
the context of superconductors. The reason frequently 
given is again the presence of the lattice, which breaks 
the translational invariance for the electrons. (Even with- 
out the lattice, a uniform background charge breaks the 
Galilean symmetry.) This is not convincing, as the trans- 
formations may well be taken to include the ions. Then 
the usual invariance under translations, rotations and 
Galilean boosts is restored, and we retain local conserva- 
tion of (total) momentum and angular momentum. So, 
as in superfluids, Vg is a perfectly legitimate thermody- 
namic variable in superconductors. It characterizes the 
spontaneously broken phase symmetry, and (as we shall 
see) gives rise to equilibrium mass and charge current 
simultaneously. 

In contrast to the exact gauge symmetry, the Galilean 
transformation property of w^, and hence Eq(|^), are ap- 
proximative and subject to relativistic corrections. This 
has been carefully considered by various authors, obtain- 
ing values 1^ that lie between 10"'' and 10~^. The rele- 
vant quantity here is the work function W of the electron, 
or the electrochemical potential, over mc^ . As pointed 
out by these authors, this value is subject to further 
scrutiny, as possibly important effects such as the cen- 
trifugal force of the rotating frame or the presence 
of the symmetry breaking lattice j^], have not been in- 
cluded. 

Eliminating these sources of uncertainties, especially 
by including the ions explicitly, the accuracy is shown in 
this paper to be given by /ichcm/c^, where /ichcm is the 
chemical potential of the superconductor, or the energy 
per unit mass of adding a (neutral) atom to the crystal. 
Due to the much greater mass of the atoms, /ichcm/c^ is 
tiny, of order 10~^°. This makes the correction negligi- 
ble for all practical purposes, and renders Eq(|l|) about 
three orders of magnitude more precise than any present 
experimental technique to determine the electron mass 
directly §. 

To understand the connection between Eq(|l|) and the 
Galilean transformation property of Vg , consider the Lon- 
don equation, which for more than half a century has 
been very useful to account for static properties of super- 
conductors, especially the MeiBner effect and the associ- 
ated magnetic healing phenomena. It is obtained by tak- 
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ing js = QsVs in the static Maxwell equation VxH = 'yjg, 
where Qs is the superfluid density, a stiffness constant. 
(The superconducting number density = Qs/m is 
probably more familiar.) Applying a further curl on 
both sides, we arrive at (what we for simplicity shall 
here call) the London equation, V x V x H = — £)s7^B, 
(B = V X A,) which yields the magnetic penetration 
depth A = ±7y^ (for H = B) ||]. Unfortunately, the 
same equation also insists on vanishing S-fields in the 
bulk, irrespective of a possible rotation ^ 0. 

Because of its curious ineptness to describe rotating su- 
perconductors, the London equation is not usually con- 
sulted in this context, (a notable exception is ^). In- 
stead, Eq(||) is considered alone, with a large dose of 
healthy intuition to compensate for its incompleteness: // 
we rotate a hulk cylinder of superconductor, it is essential 
that over the sample as a whole Vg = uiatticc? because oth- 
erwise very large currents would flow. • . Q . This appears 
convincing, as taking the curl of Eq(^) with Vs — wiatticc 
yields Eq(|^) with — ■^V x uiatticc, but it is not re- 
ally consistent: Inserting Vs = wiattico also in jg = ggVg 
leads exactly to the large electric current that was to be 
avoided in the first place. And setting Vg = wiattico while 
considering linear motions of the superconductor plays 
further havoc with physics. 

This paper reports the thermodynamic derivation of a 
generalized London equation that is (i) valid in an arbi- 
trary incrtial frame, and (ii) given in terms of the macro- 
scopic, measured fields. It is 



V X V X Ho = -gs-ri-fB/r, + 2n), 



(3) 



where Hq denotes the field in the local rest frame of the 
lattice, and 



77 = 1 -I- /ichcm/c^ 



(4) 



is the relativistic correction. Clearly, this equation ac- 
counts for the bulk value of the field, B — —2rjVl/^, and 
the attendant magnetic healing, both with and without 
rotations. 

Eq(^) is obtained by combining 



V X Ho 7js, 



(5) 
(6) 



Both will be derived carefully below, but let us first es- 
tablish a qualitative understanding: Eq(^ is an Euler 
equation, expressing the fact that the free energy is min- 
imal with respect to variations in the vector potential 5 A] 
Eq(^ is a constitutive relation, expressing the thermo- 
dynamic conjugate variable = de/dvs as a function of 
the two velocities, Vs and u„, where u„ is the velocity of 
the lattice points. As the left side of Eq(||) is invariant 
under a Galilean transformation, so must the right side 
~ ja be, which may therefore only depend on the ve- 
locity difference. Eqs(||) states this fact to linear order, 
including the relativistic correction. (This is where the 
information that is a velocity enters the equations.) 



To understand the transformation behavior of Vs, we 
start from the rest frame form of the Joscphson equation. 



(7) 



Consider first the electrochemical potential /i_ : Take the 
thermodynamic energy density e either as a function of 
the mass and charge density, 

de = ^dp + <i>dpe, (8) 

or as a function of the numbers of ions and electrons, 

de = /i+ dn+ + /i_ drt_. (9) 

These two pairs of independent variables are related by 

g^Mn^+mn^, pe = |e|(n+ — n_), (10) 

with M, m as their mass and |e|, e as their charge, respec- 
tively. As a result, the conjugate variables are related as 



M4 



/i_ = m/i + 6$. 



(11) 



Note two points: (i) This is the only place where m enters 
our considerations independently from Eq(|^), yet since 
g, Pe,n^,n- are all strictly conserved, M,m,e are the 
bare parameters, (ii) Following the standard notation 
in relativistic physics, we take the energy density e to 
include the rest energy, gc^. Then all three chemical po- 
tentials contain a large constant term, especially 

M = Mchom + c^- (12) 

With the Joscphson equation now given as 

{fi/2m)if + ^ = -e^/m, iJs + VyU = eE/m, (13) 

we may define a 4- vector 

Ua = {h/2m)da(p - -fAa = (uq, V s) , (14) 

such that a transformed u'^, in the system boosted by v, 
is given as (uq, v'^) = {uq + v • f s/c, f s + uqv/c), or with 

(|T3|,|T^ 



(15) 



Taking fj,/c^ = 1 + Mchcm/c^ = Vi see that Vs/rj trans- 
forms as a velocity, rendering the combination in Eq(|^) 
invariant. 

As is clear from Eq(^), ^chcm is the energy of adding 
an atom to the crystal at constant charge, divided by the 
mass of the atom. Estimating this energy as IC^K (boil- 
ing/melting temperature), or 10~^^J, and the atom as 
having 50 proton mass « 10^ electron mass « 10~^^kg, 
we have |/.tchem| ~ 10^(m/s)^, or 77 « 1 — 10^^". 

The main correction being so small, any inaccuracy 
m Eq(0) will be even less important. For instance, 
the velocity dependent terms in the Joscphson equation 



2 



(not considered above as the rest frame form was cho- 
sen) are smaller if the macroscopic velocity stays below 
y/JIchcm ~ lO'^m/s. The same is true for the centrifu- 
gal force, which slightly compresses the outer rim, and 
decompresses the center, of the superconductor. As a 
result, the density g, and hence /ichom(f?), arc inhomoge- 
neous. However, since /ichcm — (J^ x r)^/2 remains con- 
stant this correction is again of order (u/c)^. 

The reasons previous authors arrived at much larger 
relativistic corrections vary: Some considered the 4- 
vector Wa = {h/2m)da^ that transforms as w'^ = Ws — 
v(/x_/mc^). (See for instance Eq(lO) of 1^.) Mainly 
because of the smallness of the electron mass m, the 
associated relativistic correction /x_/mc^ deviates more 
strongly from 1. Others did consider the 4- vector of 
Eq(|lJ), especially Cabrera who in fact obtained the 
combination, /i_ — 6$ as the relevant factor in the boost 
transformation. The absence of ions in his considera- 
tions, however, prevented an interpretation of this combi- 
nation as the chemical potential of the solid. In addition, 
the term e$ was taken to be given only in the rotating, 
and not in the laboratory frame, and was therefore ne- 
glected. 

We now proceed to prove Eqs(^ by thermodynamic 
considerations, in a way that is similar to the minimiza- 
tion procedure of the Ginsburg-Landau functional JTH , 
though the procedure is executed in a general inertial 
frame here |T^, not necessarily that of the lattice. The 
basic Gibbs relation is 



de ^Tds + fidg + Vn ■ dg*°* + aijdVjUi 
+Eo-dD + Ho-dB+j, -dt;,. 



(16) 



The first four terms are that of a normal solid ||l3|], with 
^! s, p, (7*°' denoting the densities of energy, entropy, mass 
and momentum, respectively, while Ui is the displace- 
ment vector. The conjugate variables are defined by the 
respective derivatives, eg T = de/ds, or w„ = 9e/%*°*. 
The next two terms account for the presence of fields, 
where the subscript o denotes the respective field in the 
rest frame d„ — 0. The last term appears in supercon- 
ducting or superfluid phases, with the variable given by 
Eq(^. Accepting that the energy is given as in Eq(16), 
minimizing the free energy while keeping the conserved 
quantities constant, stationarity with respect to varia- 
tions in the vector potential, 

S J{e - Ts) - y"[Ho<5(V x A) - j,<5(7A)] + . . . 

quickly leads to Eq(||) . The other equilibrium conditions 

are @: 9,1;^" + djV^ = 0, dtv"- + Vn = 0, VjCTy = 0, 
Eo = 0. 

Although the basic structure of Eq(|l6|) simply states 
the fact that the thermodynamic equilibrium of a su- 
perconducting crystal depends on the specified variables, 
there are perhaps three points that may seem puzzling at 
first: (i) Why is the term ^dpe of Eq(||) substituted by 



its partially integrated form, Eq • dD of Eq(|l6[)? (ii) Why 
do the two fields, E and H, assume their local rest frame 
values? And (iii), to what extend is v„ the velocity of 
the lattice? The answer to the first question is connected 
to locality: Eq is a local function of D, while pe from 
everywhere is needed to calculate $. Local thermody- 
namics and hydrodynamics can only deal with quantities 
that preserve locality. The second question is discussed 
in details in |14|, and is a result of the fact that the 



conserved, total momentum density g contains field 
contributions. The third question is discussed in p^ , in 
which the dynamics of a hypothetical superfluid crystal 
is derived, where the equation of motion for the displace- 
ment vector Ui is shown to be iii — u" + , with 
a dissipative term that accounts for diffusion of defects 
and interstitials So, if Y^ is zero (especially true 

in equilibrium), the lattice points move with Vn- This 
remains valid also for superconductors, the dynamics of 
which reduces to that of a superfluid crystal in the limit 
e ^ of Eq(|). 

We now revisit Eq(^: When deriving the constitutive 
relation for jg , the information that it must be an invari- 
ant quantity was a consequence of asserting the consis- 
tency of Eq(^) . More prudently, this should be obtained 
as a result. So, the form of jg will be derived indepen- 
dently below, without the reference to Eq(|^), or even 
the input that Vs/rj is a velocity - though the equiva- 
lent information of Eq(|l3|), the Josephson equation, is 
needed. As we shall see, this calculation will in addi- 
tion provide further insights in, and understanding of, 
the properties of superconductors. We start the deriva- 
tion by observing that the energy current Q contains the 
term j^/^, then deduce the form of the total momentum 
density by the known symmetry of the energy stress 4- 
tensor, = Q/c^; and finally, we obtain the explicit 
form of js via a Maxwell relation linking to j^. 

The energy current Q is obtained by evaluating e via 
Eq(|l|), and requiring that it is given as a total diver- 
gence, £ = —V • Q. In the rest frame Vn — 0, and disre- 
garding dissipative terms, we have s, Ui = 0, hence only 
four terms remain: 



With Vs given as in Eq(^^, the Maxwell equations and 
the continuity equation must contain the corresponding 
counter terms. 



B = -cV X E, 
D = cV X H — ejs/m, 
e-f V-js = 0, 



(17) 
(18) 
(19) 



such that the energy current is Q = jsP -I- cE x H. The 
same calculation in a general frame t»n ^ is tedious, but 
not more difficult. And of the new terms, all Vn, the 
overwhelming one is that associated with the rest mass, 
gc^Vn- So the total momentum density is 

= gvn + i,p/c^ + E X H/c, (20) 
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while the neglected terms are again relativistic correc- 
tions, ones that are completely irrelevant in the present 
context. Now consider the Maxwell relation, 



(21) 



Evaluating the first expression via Eq(|20|) we arrive at 
the third expression; and the second equal sign only al- 
lows a functional dependence as given in Eq(|). This 
concludes the independent proof. 

The total momentum density, Eq(po|), may be written 
as — {q— Qs)vn + QsVs for E implying a literal 
nonclassical rotational inertia for superconductors - 
though it is of course much smaller than in He II: Em- 
ploying Eqs(|^,^), one finds in a superconducting cylinder 
of radius i? 3> A, that the superfluid velocity Vs devi- 
ates from the normal one, t>„ = 17 x R, by the amount 
X{2Q + jBey^t) exp[{r — R) / X], [r is the distance from the 
center, and Boxt the external field). Therefore, the frac- 
tion of Qs/ Q of the total mass does not quite participate 
in the rotation of the rest, {q— Qs) / Q, and reduces (or en- 
hances) the moment of inertia accordingly - as a function 
of the external magnetic field. 

Eqs( p^ , |l9| ) show that = de/dvg transports both 
electric charge and mass in equilibrium, with a quotient 
given by the microscopic parameter e/m. As is clear from 
the derivation, this is intimately related to the form of the 
Josephson equation, and to the nonclassical rotational 
inertia - or that 5'°' ~ js- If the Josephson equation, 
Eq(p^) , were of the form {h/2m)(p — — e$/m, the terms 
~ js in Eqs(^ 20) (ie in the momentum density and 
mass current) would vanish; if it were (h/2m)(p -f = 0, 
the persistent electric current, in Eq(|l8[), is zero. The 
first is a hypothetical case of pure superfluidity in the 
electric charge, the second a case of pure superfluidity in 
the mass, such as realized in He II. The specific sum in 
Eq(p^) characterizes a phase with condensed electrons, 
ones that carry mass and charge of the given ratio. 

All this is reminiscent of the relative broken symme- 
tries in the superfluid phases of ■^He - especially '^He- 
Ai, in which the condensation is in the up-spin popula- 
tion, of '^He particles that carry both spin and mass, with 
the ratio h/2m. And the corresponding broken symme- 
try is a linear combination of spin and phase symmetry 
- if one forgets the less obvious orbital one. The micro- 
scopic coefficient h/2m also characterizes the ratio be- 
tween the equilibrium currents for mass and spin, and it 
enters the relevant Josephson equation which - as 
in Eq(p^) - contains two thermodynamic conjugate vari- 
ables, fj, = de/dg and u: = de/ds. This gives rise to a 
number of useful analogies that should be explored. 
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